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Abstract. We analyze the first-order nonhomogeneous differen- 
tial equation -zf'(z)+Pf(z) - E(a + l;S;z)E(a + j;fi;-z) on C, where 
E(-) denotes the MacRobert's ^-function. We show that f can be 
represented in terms of the double series ^■'i'i> a ^ so known as the 
Kampe de Feriet function. As an application, we demonstrate that 
the integrals of type J P t@ +l E(-;i/t)E(-;-i/t)dt, where p e U is arbi- 
trary and 1-0,1,..., are given in terms of f . 
MSC(2010): 33C15-33C20-33C60. 



1. Introduction 

In the first part of the present letter we give some basic properties 
of the solution to the first-order differential equation, 

-4- + -) f(z) = ~E(a + l;8;z)E(a + y; /3; -z) (1) 
az z ) z 

(z,a,/3,Y,8 g C are such that the right-hand side exists), where E(-) 
denotes the MacRobert's ^-function [Mac54, Mac58l 



i 

E 



a 1 ,a 2 ,...,ap_\ n^Ha,) ^ ^ ai ,a 2 ,...,a p _V. ^ (9g) 



l&i,6 2 ,...,V *J U q ,T(bj) pFq [bi,b 2 ,...,b q ' z 

7 — 1 J 



2ni Jb 



z^nonUnaj-o 

J ~ d( (|argz|<;r) (2b) 

is n q j=1 nbj-o 

({dj g C: j = 1, ...,p}; {bj g C: j = l,...,q};z g C), where T(-) is the 
gamma function, p F q (-) is the generalized hypergeometric function. 
The Barnes-type contour B is taken up the Im^-axis to ensure that the 
poles of {T(cij-0'- j 1 = an d {T(bj-()'- j = Q) he to the right 
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of the contour and the poles of T(() to the left of the contour. In par- 
ticular, = M(-;-;-) where M(-) is the confluent hypergeometric 
function of the first kind (Kummer's function) [AS72 , §6, §13]. 
We identify f: C — ► C by a convenient notation 



a, y 

and show that it can be represented by a double series 



(a, y )_ ^ z- v T(a + j + v) 
[P, 8> Z \-L ~ " A 



v=0 



2ni Jb 



v\T(p + v) 

zenana + 1-0 



'a + l, p + v 

s, p + i + v' 



T(6-0 



E 



T(a + l)T(a + r) ^i-i i 

b 1:1;1 



'a + y, p-( 
{ P, P + l-C 

P: a + l; a + y; 
P + l: 5; P; 



1 1 

z' z 



(3a) 

(3b) 
(3c) 



Tip + DTid) 

where the integral representation d3bl ) is obtained from the series rep- 
resentation by substituting d2bl ) in 03at > and summing over v. In ( [3c1 >. 
the Kampe de Feriet function .'i.'i of two variables is defined in agree- 
ment with HGS91L 



i.l:l;l 
1:1;1 



a 1 

ci 



bi; 
di; 



d, 



Zl,Z2 



mi m<i 



(&l)mi+m2 
=0 (Cl)m 1 +m 2 (dl) mi (d2)m2 m l ] - m 2 ] - 



(4) 

i(a) m the Pochhammer symbol). If we substitute ( l2"al > in d5al > and then 
express 2F2 as an infinite series, equation d3cl > holds by ©. 

We also give the asymptotic expansion of f and some recurrence 
relations. 

Although the first-order nonhomogeneous differential equation O 
is of a rather special type, in the second part of the present letter we 
show, as an illustrative example of possible application of functions f, 
that 



rp 

I t p E(a;p;i/t)E(a + y;P;-i/t)dt 
Jo 



= p p+1 E(a;P;i/p)E 
-ip P+2 P 



'a + y, p + 1 i 
, P, P + *~~p 



a. 



7 



[P + 2, p + V p \ p 



/3+3 



a, y+1 i 
[P + 3, P + Vp 



(5) 
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where a,B,y e C; Re B > -1; p e R is arbitrary. In terms of the Rum- 
mer's function M, the integrals of type © are frequently dealt with, 
for example, in atomic physics when one considers the normalization of 
regular Coulomb wave functions [AS72, §14] of the radial Schrodinger 
equation with angular momentum I and the strength of the Coulomb 
interaction g. In this case parameters a = I + 1 + ik, B = 21 + 2 and 
y = ±2ik, where k ~ \g\ is real and positive. Another equivalent ex- 
pression of © is given in d2^1 ). 

Note that © is divergent as p — ► ±oo. Some generalization of © is 
examined as well; see ( I28D . By virtue of ( 12"B1 ) we show another proof of 
the integral ti2"§l > that has been studied in HSH031 . 

Throughout, the parameters a,B,y,8 e C and the arguments z e C 
and p,t g R are assumed to be such that the hypergeometric functions 
as well as the gamma functions are well-defined unless additional con- 
ditions are imposed. 



2. Solution 

In this paragraph we prove that d5al > solves (ffl). The equivalence 
of the remaining two expressions, d3bl-([5cl>. has already been shown 
above. 

Proceeding from ( |2aD one deduces that 



d „ 
z—E 

dz 



a + 1, B + v 
{ 5, p + l + v'~} 



. a ._,/a + l, B + v 



■E(a+l;6;z). 

(6) 



Apply zd/dz to d3al >. then substitute © in obtained expression and get 



d_ 

dz 



a > r.. 

1/3, ) 



E 



z v r(a + y + v) 



£i(v-i)ir(/j+v) 



~ 0- v r(a + y + v) d „ 
+ ) ■ z — E 

£ v\np+v) dz 

1 ( a, 



a + 1, B + v 
[5, B + l + v ; 

'a + 1, p" + v 



{ 6, B + l + v 



jS + l, o 



7+1 ^ S 2" v r(a + 7 + v) 



f vT(/3 + v) 



-S(a + l;5;z) + (^ + v)£ 



a + 1, 



B + v 



{ S, B + l + v' )) 



a, 7 



\z 



+ E 



z v r(a + y + v) 



\{v-l)\T{B + v) 
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x E 



'a + 1, p + v 
S, p + l + v 



= - E(a + 1; 6; z)E(a + y; /3; -z) + /3 



a, j 
1/3, S 



\z 



as required. 



3. Convergence 
By expanding d3ab or d3cl>-(l4l> explicitly 



T(a + l)T(a + r) + (T(a + l)T(a + r+l) T(a + 2)T(a + y) \ 1 



T(P + 2)T(8) (p + l)T(l3)n8+l)l z 



T(a + 2)T(a + r+l) 1 ^ , 

—^+U\Z ) 



(7) 



(j6 + 2)r(yS + l)r(5 + l) z 2 

we deduce that the series converges uniformly for \z\ > 1 (the parame 
ters are such that gammas are well-defined). 
Indeed, let 



(a r \ r(a + 1) 00 

r x> z \ = tv^ E »v(z), where by (HU) and (ial», 



. x z v T(a + j + v) 
Uv\z) = — r^TTT— : — - 2-F2 



'a + 1, /3 + v 1 
k 6, p + l + v'~^ 



(8) 



vir(^ + i + v) 

We show that for all e > 0, there exists N - N(e) such that, for v > N 
and for p = 1,2,..., 



\u v+ i(z)+ u v+ 2(z) + . . . + u v+p (z)\ < e for all \z\ > 1. 
One can choose N large enough so that for v> N, 



(9) 



r(v) ~ V2nv(v/e) v 
(the Stirling's formula). But 



2-^2 

for v large, and so 



'a + 1, 6 + v 1) ... . „ , 

k 5, ^i + v ; -;)~ M(a+1;5; - 1/z) 



g/3+] = y M(a + l;o;-l/2)g- v e v v g+r - /i - v - 3/2 . (10) 
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As seen, |u v (z)| < |w v (l)| and u v (z) — ► as v — ► oo for |z| > 1. By applying 
the Minkowski's inequality to © and replacing v by v + s in ( 110D for 
s = l,2,...,p, yields 



p 

s=l 



p p 

< Y, \U V+S (Z)\ < Y I"V+8<DI 
s=l s=l 



1 



( e/v )V+/3+l-a-r v -l/2 M(a + 1; 5; -1) 



f>/v) s 

s=l 



V! 



(e/v)^ 



+ l-a- 



r M(a + 1;S;-1) 



<e (v large and |z| > 1) 



for all p = 1,2,..., where e > can be chosen arbitrarily small. This 
proves that the series ( !3aD converges uniformly for |z| > 1. 

Since 2^2(0 is entire, that is, its radius of convergence is 00, one 
deduces that the series is convergent for < \z\ < 1. We note that all 
other regions of convergence due to the parameters a, ft, 7, 6 £ C can be 
obtained by applying the properties of 2-^2(0 to d3al >. 



By ©, 



a, r. 



4. Asymptotic expansion 



\ T(a + l)T(a + r) _n 



00. 



(11) 



We also note that the analysis of ([1]) yields the same result. Indeed, 
functions M(-;-;±l/z) —- 1 as \z\ — ► 00, and so by ( I2"al > the right-hand 
side of © becomes [r(a + l)r(a + yVdX/^rXtf))]^ 1 up to 0(z" 2 ). By 
solving the modified differential equation one obtains flllD . 
By Paris IPar051 . 



2^2 



a + 1, 



£ + v 



- ~(-l/^) c 



(i8 + v)r(5) 



^5, ,6 + 1 + v' j ' T(a + 1) 

for I arg z\ < n/2 and Re z < 0. Substitute ( fl2l > in © and get 



as |z|-0 (12) 



a, y 



z ~(-l/z) 



«^ 3 -i/ z g^" v r(a + r+v) 



t vT(/3 + v) 



r(a + r) (-l/z) g - <5 e- 1/2 M(a + y; j8; 1/z). 

r(j8) 1 H 



But [ASH §13] M(a + 7; # 1/z) ~ (-z) a+ iT(/3)/r(j8 - a - 7) as |z| 
Rez < 0, and so 



Ofor 
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a, y 



e -v H _ z) r +5 . Ha + r) 



T(p-a-y) 



(13) 



as |z| — >■ 0. Note that ( TTH1 > as well as the asymptotic expansion of f for 
Re z > can be established by a direct analysis of (Q]>. In the latter 
case the right-hand side of© is modified to e 1/z z^~ r T(a + l)/T(d-a- 1) 
(\z \ — • 0;Rez > 0), and thus 



a, 
ft 



r. 



1/z 0+2-y r(a+l) 

1X5 - a - 1) 



as |z|^0(Rez>0). 



5. Some recurrence relations 



By using the recurrence relations for ^-functions [in [Bho62] sub- 
stitute p = q = 2, a\ = a + 1, a,2 = ft + v, pi = 8, p2 = ft + 1 + v and get 
(fliat : in HErd53I §5.2.1, eqs. (7) and (9)] substitute a 1 = a+l,a 2 = P+v, 
Pi = 8, P2 = j6 + l + v and get ( 114bD -( fT4"cl )l one obtains from ( I3"al ) that 



£(a + l;5;z)£(a + y;/3 + l;-z) = 
P~a-1 





£-5+1 
)8 



a, 

A 


5 ;2 ; 




a, 
U 


H 


i * 



1 



z + 



a+1, y-1 
jS+l, 5+1' 

'a + 1, y 
/3 + 2, o + l 



a, y 
5-1' 



z + 



j8z s 



'a + 1, y 
^ + 2, o + l' 



(14a) 
;z) (14b) 
(14c) 



Since the proofs are obvious due to the properties of ^-functions we 
will omit them by noting that the recurrence relations d!4bl )-( fT4c1 ) are 
reduced by applying dl4aD so that the sum 



~ z- v Y(a + j + v) E 
v to v\T(p + v) 



'a + 2, p + l + v 
5 + 1, /3 + 2 + v 



\z 



can be calculated by setting r(/3 + v) = T(/3 + 1 + v)/(/3 + v), thus yielding 

a+1, y 



(a + 1, r -i i 



,;z 



^ + 2, 5+1 

where the first summand is then represented by d!4aD 
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6. Some series for ^-function 
Here we prove the sums arising in the subsequent paragraphs: 
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and 



oo (_iy 

s ti(j0+D s 



(a + l, [2] s _i 
j3 + l, [l] s -i' 



z\ = -E 



(-l) s (r + s)! 



-(P + l) r+1 E 



E 



a + l, [2] s _i 



1)3+1, [l] s _! 

fa, [jfl + lUi 
1/3, [j8 + 2] r+ i' 



a+l, /3 + 2 
1)6 + 1, yS + 3 ; " 
i 

zr\ E(a;P;z) 
(r = 0,l,...), 



S ( - i),£ (^: 



a + l, [2] s _i 



[1] 



s-l 



;z\E 



-P 



U5 + 2, /3+r z J 2 



a + y, [j8 + l] s+ i 
k ft t)6 + 2] s+ i' 

7+ 1. 



a, 



1)3 + 3, J8+1' 



(15a) 



(15b) 



(15c) 



where [a]„ (a e C;n = 0, 1,...) denotes, for short, the set of parameters 
a, a,..., a (n times). 

Apply d2bl > to the left-hand side of ( I15al > and write 



OO (_1)S 



'a + l, [2] s _i 
,13 + 1, [1] 



l' / 2ni Je 



zenana +i-o 



x °° (-i)'(i-O'- 1 
x s ti (p + iy 



2ni Jb 



b np+i-o 

z<nQT(a + 1-0 1 



np+i-o p+2-( 



By applying <|2b) along with (/3 + 2 - = T(/3 + 3 - OWP + 2 - to the 
latter contour integral yields d!5aD . 

Similarly, apply d2bl > to the left-hand side of d!5bl > and get 



™ (-l) s (r + s)l E 

s ti sKp+iy 



'a + l, [2] 8 _i 
J0+1, [l].-i' 



27ri Jb 



/rxprxa + i-p 
r(js+i-o 



xE 

s=l 



S (-l) s (r + s)!(l-O s_1 1 f 7y 2 f r(0r(a + l-0 r! 



-f 

2ni Jb 



B np+i-o w 



-i+ 



s!(/3 + l) s 
(P+D r+1 

(p + 2-O r+1 ) 2ni J B ^ r(j6 + l-0(l-0 2^i 

2^r«)r(a+i-o 



_r!_ r /r(or(a + i-o + r!(/3 + ir +1 

2777 J/J 



b T()6 + l-C)(l-0()6 + 2-0 



(16) 
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In order to calculate the obtained contour integrals, two possibili- 
ties are valid: either apply d2bl and then reduce the order of ob- 
tained E'-functions or make use of the residues directly recalling that 
Res^ = _v T(0 = (-l) v /v! (v = 0, 1, . . .)• The combination of both strategies 
yields 



1 f z ( T(QT(a + 1-Q =E 
2ni Jb 

7ii Jb 



(a+1, 1 



\z 



T(a) 



■E(a;p;z)\, 



2ni jb 

oo 

E 



r(j8+i-od-0 \P +1 > 2 

/rXprXa+l-Q _ F (a + l, [fi + 2\ r+ i, 1 ] 

r(^ + l-0(l-0(j8 + 2-O r+1 [P +1 > [/5 + 31r+i, 2' Z \ 



(-l/^) v r(a + l + v) 



zT(a) 



J (v + l)!r(/3 + l + v)(/3 + 2 + v) 



r+l 



r+l 



(-l/z) v-1 r(a + v) 



zT(a) (a, [/3 + l] r+1 ) 

r(/3)(j8+iy +1 [h [j8 + 2Ui' Z J 



v t- v!r(j6 + v)(j6 + l + vr +1 

Substitute obtained expressions in ( fTBl > and get d!5bl ). 

The application of d2bl > to the left-hand side of ( |15cl > gives the inte- 
grals 



s=l 



a + 1, [2] s _i 
j8 + l, [l].-i' 



2 £ 



fa + y, [/3 + 



27ri Jb 



/r(or(a + i-o i 



E 



r(j8+i-o 
(-im-0*- 1 



[j8 + 2] s+ i' 

,(-zfnona+y-o 

CM, 



2niJ B T(p-0 



— [ 

2ni Jb 



/r(or(a + i-o 



i (j8 + 1 - C0 S+1 2ni Jb H/3 + 1 - 

(-«)<'r(Or(a+7-0 



2ni Jb 



r(/8-c , X/s+i-c / Xj8+2-c-C0 



/r(or(a + i-o 



a + r, /3 + 1, yS + 2-C 
A j6 + 2, j6 + 3-f 



r(js+i-o 

By the recurrence relation for ^-function [in [Erd53, §5.2.1, eq. (9)] 
substitute ai = a + y, a<i = P + 1, 03 = (5 + 2 - (, p\ = p + 1, p2 = P + 2, 

p 8 =/3+3-a 



(a + r, p + 1, P+2-( _ 
{ 0, p + 2, P + 3-C 

a + 7+1, /3 + 3-C. 



PE 



'a + r, P + 2-( 
+ P + 3-C 



1 

+ -E 

z 



P + 3, P + 4-C 
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and so 



£ 1 ( - 1)J V1, [1] 

2ni Jb 



np+i-o 



2n 



i r jy /r(or(g + i-o p 

— d( — — : — E 

tiz Jb 



'a + y, [/3 + lL+i _ 
k P, W + 2] s+1 ' Z 

'a + y, p+2-( _ 
j8 + 2, P + 3-C ' 

'a + y+1, P + 3-( 
k /3 + 3, p + 4-( 



np+i-o 

But these two integrals are nothing more than the integral representa- 
tions of function f, where the parameters ft, y, 8 in ( 1351 ) are accordingly 
replaced by p + 2, y, p + 1 in the first integral and by ft + 3, 7 + 1, /3 + 1 
in the second one. Therefore, equation d!5cD holds. 



7. Definite integral 

In order to calculate the integral © we make use of the series ex- 
pansion ( l2"al ) of ^-functions (equivalently apply (I3"cl)-(H1)), 



I 



p n , , , S i^(-t) v r(a + /i)r(a + r + v) 

t (i E(a;P;i/t)E(a + y;P;-i/t)dt= £ ' 



~ Mivir(j8+/i)r(j8+v) 



Jo 



(17) 



The integral on the right-hand side equals p^ + ^ +v+1 /(p + p + v+1) pro- 
vided Re/3 > -1. Expand (/3 + /i + v + l) -1 into the binomial series, 



P+p+v+1 p 



_ 1 * (-lY +s p r v s (r + s)\ 

1 ~J+i r h_ r\ s \(p + iy +s 



(18) 



and get 



I 



„/3+l 00 ( — Y) r+S ( r + c)! 

g (ipYp r T{a + p) g (-ip)VT(a + y + v) 



^-J, /i!r(j8 + M) v= vT(/3 + v) 
In fn?l ) the latter two sums are of the type 



(19) 



°° z v v r T(a + v) t A , x 
L rFT^T x (r = 0,l,...). 
^ vT(/3 + v) 
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For r = 0, this is a usual hypergeometric series given by E(a; P;-l/z). 
For r = 1, make a substitution v — ► v + 1 and get zE(a + 1; /3 + l;-l/z). 
For r > 2, write 



S 2 v vT(a + v) _ S z v (l + v) r_1 r(a + 1 + v) 



°° 2 v r(a + l + v)[r(2 + v)f- 1 _ 

: v tbv!r( ) 6+i+v)[r(i+v)F- 1 ~ zE 



(a + 1, [2] r _i 1 
j8 + l, [1W z 



(20) 



Note that ( f20b also holds for r = 1, provided [a]o is the empty set. 

In dl9D expand both sums over r > and s > into two terms, one 
with r = and s = and another with sums over r > 1 and s > 1. Then 
substitute the right-hand side of ( 120D in obtained expression for r > 1 
and s > 1, and E(a; p"; Up) and E(a + j;P;-i/p) accordingly for r = and 
s = 0. The result reads 



+i 



CP 

I t fi E(a;P;i/t)E(a + r,P;-i/t)dt = 
Jo p + 1 



£(a; /5; i/p)E(a + y; P;-i/p) 



fa + l, [2] s _i n 2 ^ (-l) r fa + l, [2] r _i H 
1/3+1, [l] s -i' p J +P r tlr!(j8 + ir 1)8 + 1, [lWpj 



oo t-iy 

s ti(i6+i) s 

°° (-l) s (r + s)! /a + y+l, [2],_i i 
x 2_ — — — E 



s=l 



s!(p>l) s 



i )8 + i, [1W p , 

The sums over s > 1 are calculated in A15al >-( fl5b~l >. By substituting 
them in the above expression, after some elementary simplifications 
one obtains the integral 



rp P +1 

/ t fi E(a;P;i/t)E(a + r,P;-i/t)dt = - 

Jo p+1 



E(a;fi;i/p) E(a + j;P;-i/p) 



+ 'H + 1, ^ + 3 ; -p)J + ^ + 1) s ? 1 ( - 1)£; 



a + l, [2] s _i r 
jS+1, [l] s _i' P; 



a + y, [/3+l] s+ i __i 
i P, [j6 + 2] 8+1 ' o 



The latter sum is given in ( !15cD . By substituting the right-hand side 
of dl5cl > in the above expression and noting that [in [Bho62> eq. (4.4)] 
substitute p = q = 2, a\ = a + y, a<i = P + 1, pi = P, pi = P + 2, and after 
that apply [Erd53, §5.2.1, eq. (7)] with ai = a + y- 1, pi = p"- 1] 
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E(a + j;(5;-i/p) + ipE 
= (P + 1)E 



'a + 7+1, (5 + 2 i 
k P + l, i6 + 3 ; "p 

'a + y, (5 + 1 i 
, P, i6 + 2 ; "p 

one derives ©, as required. 



(21) 



8. Some generalization 

The procedure to calculate the integral © suggests a natural gener- 
alization, 



rP 

\ t p+l E{a;P;ilt)E{a + r,P;-Ut)dt 
Jo 



(22) 



for some I = 0, 1, . . . In this case, the integral on the right-hand side 
of ( fT7l ) is calculated by replacing (5 by (5 + 1, for Re (5 + 1 > -1. The 
subsequent denominator in ( [TBI is replaced by 



1 



(-l) r+s v s (r + s)l 



(5 + l + p + v+l (5+l r % Q H«!(j8 + ir „ =0 
and so ( flT?l > becomes equal to 



e k r_ "' 



tP +l E(a;P;i/t)E(a + T ,l3;-i/t)dt = Z—- £ , ' 1V ' 
Jo p + 1 r!s!(/3 + l) 

x f E f r V~* E ap)At/J " r(g + /j) | £ (-^) VySr (q + r+v) (23) 



u=o 



r pirxp^ + p) v to vin/3 + v) 



For s = 0, the sum over v equals E(a+j;/5;-i/p), whereas for s = 1,2,..., 
it is given by ( 1201 ) with a and z replaced by a + y and - ip, respectively 
(r = s). 

For r = 0, the sum in the parentheses in d2Ul ) equals i£(a; /3; f/p), and 
for r = 1, 2, . . ., it is given by (after the summation over p by ( 127)1 



l r E(a;P;ilp)+ipY J \ \l r ~ u E 



u=l 



U | 



a+1, [2] u _i i . ,. 
J3 + 1, [l] B _i ; pJ =Z E{a >P> l, P ) 



. S (ip) v r(a + l + v) ' 

+ 'pE ,„, imfll r^ E 



v f (v+i)!r(^ + i + v) u tiU 



(v+ifr-^r^a^z/p) 



oo 



j P rE 



(ip) v r(a + l + v) 



(ip) v r(a + l + v)(Z + l + v) r 



v f (v+l)!r(/3 + l + v) r v to (v+l)!r(j8 + l + v) 
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= l r E(a;P;i/p) + V 
= E 



a ■, .\ £ (ip) v r(a + v)(l + vY 
—— -E(a;P; i/p)\ + £ — , r ,„ : 

r(^) ) ~! v\r(p + v) 



'a, [l + l] r i 

fi, Mr '~p, 

Combining the results all together and substituting them in d23l > we 
find that 



I t p+l E(a; P; i/t)E(a + 7; P;-i/t)dt = 

Jo P + 1 

00 

xE(a + y;P;-i/p)-ipE(a;P;i/p)Y, (R ^,, S E 
00 f-iy 

+E(a + r ,P;-i/p)ET7r-7r E 
s =i \P + *-) 



r 



E(a;P;i/p) 

'a + y+1, [2] s _i _r 
k P + l, [l] s _i 5 Pt 

fa, [l + l] s f, (-1Y 

[P, il], 'p) iP kirl(P + D r 



(a, [I + l] r i \ g (-l) s (r + g)l g (a + 7 + 1, [2] 8 _i.__i 



XjE lA [Z] r ' p j s ^ s!(j6 + l) s "I P + l, p, 

Substitute the sums 015 aft — ( fl5bl > in the above expression and get by 
applying (ED, 



/ t p+l E(a;P;i/t)E(a + y;P;-i/t)dt 
Jo 



= P 



E(a;P;i/p)E 



a + j, P + l i 
I A P + V~~p 



+ J^i-lfE 

s=l 



fa, [l + l] s i 



'a + j, [j8 + l] 8 +i _i 
k ft [/3 + 2] s+ i 5 p 

By J2bl >. the latter sum can be represented by 



(P, ui s s ''p ]E 



(24a) 



8 p 



s=i yp> 

2ni Jb 



'a + y, [j0 + l] s +i _i_ 
k p\ [jS + 2],+i' p 



(i/p/ r(or(a -0 1 r d , M/p/'noiXa + 7 - C) 

27rj Jb 



r()6-o 



r(j8-o 



(-l) s (/-Q s 

l\s+l 



i(j8 + l-f> 



-f. 



(i/p)<r(or(a-o j_ r d , 

2^i Jb 



r(j8-o 



(-i/p/ , r(Qr(a+7-Q 
r(/3-o 



+ 



p + l-C p + l + l-(-(' 
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\ P, P + z p] 2nijB 



xE 



(a + Y, p + 1 i\ 1 f ^(ilp) ( T(Ona-0 
{ P, 

(a + Y, p+l+l-( i 



np-o 



v P, p + l + 2-C p 
and so d24al > becomes equal to 



rP 

\ tP +l E(a;P;i/t)E(a + Y;P;-i/t)dt 
Jo 



p M + i , (i/ p )(T(OT(a-0 



>-t-t-t-_L n 

ini Je 



2ni J B Tip-0 



'a + Y, P + l + l-( i 
i P, P + l + 2-('>~-p- 



(24b) 



The integral representation d24bb of d22D is more convenient than d24aD 
because in the latter case one can directly apply the recurrence rela- 
tion [in [Erd53, §5.2.1, eq. (9)] substitute a\ = a + y, a 2 = P + I + 1 - (, 
p\ = P + l, p2 = P + l + 2-( and x = -i/p] 



E { p, p + l+2-C~~ l ~ PE 



'a + Y, P + l + l-( i 
j3 + l, P + l + 2-(~~ p 

a + y+l, p + l + 2-( i 
P + 2, p + l + 3-( ; ~~p 



ipE 



- (25) 



I + 1 times for a given I = 0,1,..., and then apply d3bl > to d24bl ). 

For example, if / = 0, then d24bl > coincides with ©. Namely, sub- 
stitute / = in ( l2~5l >. and the obtained expression in d24bl >. and get by 
applying W, 



rP 

I t p E(a;P;i/t)E(a + y;P;-i/t)dt 
Jo 



P 



a-1, 7+1. i 



[P + l, P 



-\-ip 
P 



a-1, 7 + 2 i 
.P + 2, P > 



(26) 



Apply d!4al > to both f functions in the above expression and get ©. 

As an example, we also calculate the integral 022D for 1 = 1. In this 
case apply d25l ) two times, 



E 



(a + Y, P + 2-(._i 
{ p, P + 3-C p' ; 



(P + DE 



a + Y, P + 2-( i 
[P + 2, /3 + 3-C~p 
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-ipE 
-ipE 



a + y+1, 


p + s 




P + S, 


P + 4 


-c 


a + y + 2, 


P + 4 


-( 


P + 4, 


P + 5 


-r 



(P + 2)E 



'a + y+1, j8 + 3-C_r 
t /3 + 3, i6 + 4-C ; p, 



Substitute the right-hand side of obtained expression in d24bl ) and get 
by d3B, 



/' 

Jo 



t^ +1 E(a; P; i/t)E(a + 7; P;-i/t)dt = p p+z /3(/3 + 1) 



/3+2 



a 



1)6 + 2, 



-2ip(j8 + l) 



a-1, y + 2_ i 



1/8 + 3, 

For arbitrary I = 0, 1, . . ., we find by 



a-1, 7 + 3 i 
J6 + 4, j8 > 



7+1 r 



a+7, p+i + i-( i\_ l ^ (-iprr(p+i + i) (i + i) 

k p\ P+1+2-? pj-j^ n/S + n) n j 



a+7+re, p+l +n+l-( i 
fi + l + n + 1, p + l + n + 2-('~~p~ 



(27) 



By substituting the right-hand side of ( I27D in d24bl ) and making use of 
we finally derive the expression 



CP 

\ t^ l E{a\P\ilt)E{a + j\P\-Ht)dt 
Jo 



^ +l+1 l j} (-i P rnp+i + i) 



T(P + n) 



i + i\< 

n 



a-1, j + n + 1 i] 
[P + l + n + 1, p >J KZ0) 



(Re P + l > -l;l = 0,1,...). We note that, if on the left-hand side of ( 
we replace (i/t) and (-i/t) accordingly by (hit) and (-A/it) for A e C, then 
on the right-hand side one needs to replace (-ip) in the numerator and 
(i/p) in f by (p/A) and {XI p), respectively We shall omit the proof since 
it is obvious. [Hint: starting from ( l2"cfl >, one needs to replace (ip) and 
(-ip) by (-p/A) and (p/A), and then proceed similarly as in the case for 
A = i.] 

In HSH03H the authors have shown that [in [SH03 , Lemma 1] substi- 
tute a = a, b = b' = p, a' = a + j, d = p + 1, k = -k' = i and make use of 
( l2"al > of the present letter] 
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H tPe- ht E(a; P; i/t)E(a + y; ft-i/t)dt = ^ T{a + ^ 

xF 2 (P + l;a,a + r,P,P;i/h,-i/h) (Re/3 > -1; \h\ > 2) (29) 

where i^O) denotes the Appell's hypergeometric function. We shall 
demonstrate that d29D is easy to obtain from H28D . 

Recalling that exp(-ht) = l^ =Q (-ht) 1 we can multiply C§9> by 
(-h) l /l\ and then perform the summation over I = 0,1,... This gives 
the integral 



/ ^e~ ht E{a; P; i/t)E(a + y; P;-i/t)dt = p? +1 Y ■ 
Jo l=o 



i-ph) 1 



11 



{-ip)^T(p + l + l)\l + l\ 

x L 



n to np+n) 



n 



a — 1, j + n + 1 i 
fi + l + n + 1, p > 



(30) 



(Re P > -1). Applying d3bl > to the function f on the right-hand side of 
( 13"0~1 ), one finds by virtue of ( 12~T1 >, 



/' 

Jo 



t p e~ ht E(a; P; i/t)E(a + y; P;-i/t)dt - P 



3+1 



f d( (i/p) ( T(QT(a - Q g (-ph) 1 E 
Jb IY/i rx 



j B np-o /=0 

By applying d2bl > to ^-function, 



2ni 

<a + y, p + l + l-( _V 
P, P + 1 + 2-C p> 



rP P +1 

I t^e~ ht E(a; P; i/t)E(a + y; P;-ilt)dt = 

Jo 2th 



(i/p^T(QT(a-0 J_ r (-i/pfT(C)T(a + r-C) 
np-0 '2niJ B T(p-C) 



*ltP!* 1 . (31) 

The sum over/ equals (phy( i - 1+(+i: '(T(p + l-(-(')-T(p+l-(-C,ph)), 
where the incomplete gamma function T(a, t) = u a ~ x e~ u du. 

Starting from now, we can either (i) take the limit p —* oo so that 
T(-,ph) —> or (ii) apply the above given integral representation of the 
incomplete gamma function for arbitrary peR. Both cases yield iden- 
tical results. 

In case (i), the integral < T5TT > becomes equal to 
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h-P- 1 f ,AihfT(OT(a-0 



tPe- ht E(a;l3;i/t)E(a + r,P;-i/t)dt = —— dt;- r/ „ ^ 
Jo 2ni Jb U/3-0 

2ni 



r , (-ihfnona + y - onp + 1 - c - o 
) B "^ np-o 



The latter contour integral is calculated by applying the residue theo- 
rem, 



— f 

2ni Je 



(-ihfnona + r - onp + 
'b np-o 



_ °° (-t//t) v r(q + y + v)r(jg + l-C + v) _ r(a + r )r(/3 + i-o 

" ho vinp+v) " npi 

x afi(a + y,j8 + l-C;j8;-i/«, 
and thus 

tPe- ht E(a; p; i/t)E(a + y; )6; = - ■ — d( 

Jo hP +1 T(p) 2tiiJb 

x (ih) ( T(OT(a - 0(P ~ aFi(o + y, + 1 - f ; j8; 
or by the residue theorem, 



/ t p e~ ht E(a; p; i/t)E(a + y; P;-i/t)dt 
Jo 



T(a + r) 00 (ilhY 
= h p + i Y{p) L — r <« + + v) 5^1(0 + y, j8 + 1 + v; 0; -i/A) 

^r(a)r(a + y) ~ (i/W v (a)v(^ + Dv „ , « - « 

= - — s-; — > 2Fi(a + r,P + l + v;P;-i/h). 

hP +1 T(p) to v ^y 

The latter sum is nothing more than the Appell's double series given 
on the right-hand side of d2^]>. The series i^O) is absolutely convergent 
for | A | > 2. Hence W holds. 

In case (ii), we obtain from d3Tl >. by applying the residue theorem as 
above, 

pT(a)T(a + j) 



np) 



[ tP e - ht E(a; p; i/t)E(a + y; P;-i/t)dt = - „\ 
Jo hP +1 

v<a i ««•/*, f (ih)^T(()T(a — 
xF 2 (p + l;a,a + r;P,p;i/h,-i/h) — — — d( — 

2ni Jb Up-O 
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1 f . rl (-ih)(T(C)T(a + r-(')T(l5 + l-(-C,ph) 

x / ac , (32) 

2niJ B T(P-C) 

where the contour integral involving the incomplete gamma function 
is calculated by using the integral representation of T(-,ph) along with 
the residue theorem, namely, 



(- ihfnOKa + r - on/3 + 1-C-C, ph) 

CLi, 



= S (- i/h) v T(a + y + v)T(P + l-( + v,ph) = f°° dtt M e -t 
v ^ v\T(p + v) Jp h 

S (-it/h) v T(a + r + v) r°° a r „ ., , x , 

x E 7 = t^e-'Eia + 7; P;-ih/t)dt 

v\T(p + v) Jph 

= h p+1 ~ ( / t p ~ ( e~ ht E(a + r ,P;-i/t)dt, 
J P 

where in the last step the substitution t—*ht has been initiated. 
Substitute (HD in dHU> and get by W, 



(33) 



th- ht E(a; p; i/t)E(a + y; ft - Wtff = 
o hP +L T{(5) 

x F 2 (P + l;a,a + y; ft, P;i/h,-i/h) 

- / t p e~ ht E(a;P;ilt)E(a + y;P;-ilt)dt. 
Jo 



'P 

But J P + J p °° = J °° for all peU, and hence d29]» holds again. 
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